Exchange effects on electron transport through single-electron spin-valve transistors 
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We study electron transport through single-electron spin- valve transistors in the presence of non- 
local exchange between the ferromagnetic leads and the central normal-metal island. The Coulomb 
interaction is described with the "orthodox model" for Coulomb blockade and we allow for non- 
collinear lead magnetization directions. Two distinct exchange mechanisms that have been dis- 
cussed in the literature are shown to be of comparable strength and are taken into account on equal 
footing. We present results for the linear conductance as a function of gate voltage and magnetic 
configuration, and discuss the response of the system to applied magnetic fields. 



PACS numbers: 85.75.-d, 72.25.Mk, 73.23.Hk 



I. INTRODUCTION 



Downscaling magnetoelectronic devices to the nanometer regime implies that electron-electron interaction effects 
become prominent, as has been amply demonstrated by many experimental studies on the Coulomb blockade in double 
tunnel junctions with ferromagnetic elementsiSiiiiSiSiLSiii Much of the theoretical work focusses on F|N|F spin valves, 
in which the island is a normal metal (N) and the contacts are ferromagnets (F) with variable magnetization directions. 
Initially, the interest was mainly focussed on the giant magnetoresistance, i. e. the difference in the transport properties 
for parallel or antiparallel magnetizationsii2iii*i2ii&ii More recently, the interplay between spin and interaction effects 
for non-collinear magnetization configurations has attracted quite some interest ^igiiLi2ii2i22i2ii22i22i2ii2Si 

A single-electron spin-valve transistor (SV-SET) is an F|N|F spin-valve with a sufficiently small normal-metal (N) 
island that is coupled to the ferromagnetic leads by tunnel barriers. When the electrostatic charging energy of the 
island is larger than the thermal energy, charge transport can be controlled at the level of single electron charges by 
varying bias and gate voltage, as is well known for non-magnetic SET's^ With spin-dependent electron tunneling 
rates and sufficiently long spin-decay lifetimes, a spin accumulation (or non-equilibrium magnetization) that strongly 
affects electron transport may build up in the nonmagnetic island. 

In this article, we discuss the transport characteristics of metallic SV-SETs in the Coulomb blockade regime, allowing 
OO 1 for arbitrary, noncollinear magnetization directions. In particular, we examine the influence of exchange effects through 
F|N tunnel contacts on the spin accumulation in the center island, presenting a more complete discussion compared 
to that in Ref. We argue that two separate exchange effects have to be taken into account. On one hand, 

there is the non-local interface exchange, let us call it "XI" in the following. In scattering theory for non-interacting 
systems it is described by the imaginary part of the spin-mixing conductance^ while in the context of current-induced 
magnetization dynamics XI acts as an "effective field" Such an effective field has been found experimentally to 
strongly affect the transport dynamics in spin valves with MgO tunnel junctions^ This effect has recently also been 
involved to explain magnetoresistance effects in carbon nanotube spin valves^ and called spin-dependent interface 
phase shifts^ The second exchange term ("X2") is an interaction-dependent exchange effect due to virtual tunneling 
O ■ processes that is absent in non-interacting systems. It has been considered for islands in the electric quantum limit, 
' in which transport is carried by a single quantized level onlyiii The X2 effect is potentially attractive for quantum 
. information processing, since it allows to switch on and off effective magnetic fields in arbitrary directions just by a 
gate electric potential. We compute here X2 for a metallic island in which size quantization is not important. We 
find that both exchange effects are of comparable magnitude and affect the transport properties in a characteristic 
■ way, but can be separated in principle by employing the gate dependence of X2. 

The paper is organized as follows. In Sec. [H]we introduce the model system for the SV-SET. In Sec. IIIII the two 
relevant types of exchange processes are discussed. Charge and spin transfer rates are determined in Sec. IIVI Finally, 
we present results for the transport characteristics as a function of magnetic configuration, gate voltage and applied 
magnetic field in Sec. [V] 
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II. MODEL SYSTEM 



An SV-SET (see Fig. [IJ,) is composed of a small metallic cluster in contact with two large ferromagnetic electron 
reservoirs in thermal equilibrium characterized by magnetization directions rti\ and with m a = (sw.6 a , 0, cos9 a ) 
(for a = 1, 2), where B x = 9/2 and 9 2 = -6/2. 
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FIG. 1: (a) The spin-valve single-electron transistor: A small normal-metal island tunnel-coupled to two large ferromagnetic 
leads. The unpaired spin angular momentum on the island is denoted by s. (b) The magnetization directions in the leads 
define an angle 8. 



The F a |N contacts are tunneling barriers with conductances that depend on the electron spin, G£J for the majority 
and GjW for the minority spin in the ferromagnet. The total conductance for contact a is then given by G a = 
(G£t + G l J) and the contact polarization is defined as P a = (Gj, T - G l a l ) / (G^ + G l a l ). The resistances R a = 1/G a 
are taken to be much larger than the resistance quantum Rq = h/e 2 , which, at low enough temperatures and bias 
voltages, allows us to study the blockade of transport by the Coulomb interaction. The electron tunneling rates are 
governed by the change of electrostatic energy of the whole circuit upon transfer of an electron. The capacitances of 
the junctions C a determine the charging energy of the island. 

We limit our considerations to islands formed by metallic clusters for which the thermal energy (fcsT) is much larger 
than the average single-particle energy separation (reciprocal density of states) S = 1/pn, but much smaller than the 
single-electron charging energy. Therefore, many levels on the island participate in the transport and we may treat 
the electronic spectrum as continuous. For a gold cluster with a diameter of 10 nm, S approximately corresponds to 
a temperature of 2 K. The Kondo physics of quantum dots connected to ferromagnetic leads^L22iS is suppressed in 
this regime. 

Since the currents flowing into and out of the cluster are spin-polarized, the island may become magnetized. The 
number of unpaired spins on the island is limited by spin-flip scattering, which we parametrize by a spin-flip relaxation 
time r s f. There is evidence from several experiments that the spin-flip times in metallic nanoparticles can be much 
longer than in bulk systems, which implies that the effects of a spin accumulation on the island should be taken into 
accountiS&S For later convenience we introduce the spin-flip conductance parameter G s / = p^e 2 / (2r s f). We assume 
that the energy relaxation on the island is fast compared to the electron dwell time. 

The total Hamiltonian for the SV-SET is 

H = H N + {H Fa + H Ta + H cxa ) , (1) 

a=l,2 

where Hjy is the Hamiltonian for the normal metal cluster in the "orthodox model" for Coulomb blockade 

tt t , e 2 (n N - C G V G /e) 2 

Hn = 2^ £ kc\ s c ks + — . (2) 

ks 

Here ct is a creation operator for an electron state with orbital index k and spin s € {|, J,}, where the z-axis is chosen 
as spin quantization axis. The Hamiltonian includes an electrostatic interaction energy which depends on the junction 
capacitances C a , the gate voltage Vq and the excess number of electrons on the island um- The gate voltage shifts 
the potential and induces a charge CgVg- We assume that the gate capacitance Cq <§C Gi,G2, and in the following 
we use C\ = C2 = C/2. The energy levels in the two ferromagnetic leads (denoted by a = 1,2) are spin-dependent: 

Hpg = ^^^aksO'aksO'aks- (3) 

The operators a^ fcs create electrons with spin s in the spin-quantization axis along m a . 
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It is convenient to introduce annihilation operators c a k s for electrons in the normal metal defined for a quantization 
axis in the direction of rrt a . The relation between operators in the two bases is then c a k s — U ss i (8 a ) C}. s i, expressed 
in terms of the spin ^ rotation matrix 

U (9 a ) = = ( COS y2 sin£ > Q /2 \ 

y ' V — sm6a/2 cos# Q /2 ) v ' 

Then, for each contact, a tunneling Hamiltonian 

H Ta = J2Ks a ik s ^ qs + h.c. (5) 

kq.s 

describes the coupling to the island. The tunneling coefficients are assumed to not significantly depend on energy on 
the scale of the charging energy. We discuss the exchange contribution, represented by the Hamiltonian H cxai in the 
next section. 

III. EXCHANGE EFFECTS 

Here we discuss two different exchange effects that affect the electrons in the normal metal island attached to 
magnetic contacts. These two flavors arise when the transport properties for an SV-SET are determined to lowest 
order in the tunnelling probabilities. 

A. Nonlocal interface exchange (XI) 

The non-local exchange coupling between ferromagnetic films through a normal metal spacer is an important effect 
that determines the ground state of magnetic multilayers (see Ref. |3J for a recent review) . Electrons in a normal 
metal that are reflected at a contact to a ferromagnet, pick up a phase depending on the electron spin relative to the 
magnetization direction. In sufficiently clean and narrow F|N|F structures, quantum well states are formed in N whose 
energy depend on the magnetic configuration through the spin-dependent phase. By a rotation of the magnetization 
directions the energy spectrum and Fermi energy varies, causing the ground state energy to depend on the relative 
angle 8. In metallic multilayers with a suitable spacer thickness, this can lead to an antiparallel ground state, 
which displays the celebrated giant magnetoresistance when the magnetizations are forced into a parallel direction 
by a magnetic field. Even when the ground state energies are not significantly affected by the exchange coupling, 
configuration-dependent quantized states can still be observed in transport. This has been shown for high-quality 
planar tunnel junctions^ as well as spin valves in which the node is formed by single carbon nanotubes with a 
quantized energy spectrumii^& In Ref. the effect of interfacial phase shifts on the magnetoresistance of ballistic 
quantum wires between ferromagnetic leads was calculated. The spin-dependent phase shifts give rise to a slightly 
different quantization condition, which can spin-split the energy levels. Since we are here interested in classical islands 
with a continuous electron spectrum, we calculate energy shifts for a semiclassical island using the Bohr-Sommerfeld 
quantization rule in Appendix IaI 

Here we consider the limit of tunnel junctions between a normal metallic island and ferromagnetic electrodes. The 
torques on the ferromagnets are then very small. The exchange coupling does not significantly disturb the ferromagnets 
in this limit, but persists to affect transport. The present study focusses on the charge transport properties in the 
limit of small tunneling matrix elements, thus from the outset excluding resonant tunneling, co-tunneling or Kondo 
type physics. The states on the island may be size quantized, i.e. the energy level spacing exceeds the thermal 
energy ("quantum dot"), or, in the opposite limit, better described by a semicontinuous density of states ("classical 
dot"). Here we concentrate on the latter, i.e. semiclassical, diffuse, or chaotic islands, for which it can be shown 
quite generally that equilibrium spin currents are suppressed^ The state of the island is then characterized by a 
semiclassical charge and spin distribution function that has to be determined self-consistently as a function of the 
junction parameters and the applied voltages. For non-interacting systems, the spin and charge currents through an 
F|N interface are determined not only by the conventional conductances G£J and G+* introduced above, but also by 
the complex spin-mixing conductances G^-jSl which are discussed in Sec. IIVI The real part ReG^ is the material 
parameter that is proportional to the spin-transfer torque^2i^ The imaginary part Im G]} reflects the spin-dependent 
interface phase shifts and affects the magnetization and spin accumulation dynamics as an effective exchange magnetic 
field parallel to the magnetization directioni^2i2S4£ Im G]} is relatively small for intermetallic interfaces^ but is in 
general comparable in magnitude to the other conductance parameters^S*^ The non-local interface exchange has 
been discussed in similar terms for spin valves consisting of Luttinger liquids with ferromagnetic contacts^ 
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The blocking of transport by the Coulomb charging is usually described by Fermi's Golden Rule (see below), which 
employs a probability (squared matrix elements) and energy conservation. As long as the charging energy is much 
smaller than atomic energy scales (like the Fermi energy), the junction parameters such as the interface transparency 
and spin-mixing conductance are unaffected and the Coulomb blockade is governed by the energy conservation criterion 
only. This implies that the exchange effect can be described by the Im G]} of the bare junction. 

It remains to parametrize the exchange in the limit of the tunneling Hamiltonian, i. e. to lowest order in the interface 
transmission. We show below that this is achieved by adding the following exchange term H cxa to the Hamiltonian 
for the two leads: 

^cxq — ^ ^£aksC a fc s Caks • (6) 
ks 

The energy shifts Ae a k s , see Eq. I|A2I) . are proportional to the inverse density of states , but they remain relevant for 
small level splitting because the dwell time is inversely proportional to the average energy level separation or inverse 
density of states 5 — p^ 1 . This Hamiltonian is an effective Zeeman splitting caused by an exchange magnetic field in 
the direction of the magnetization, see Sec. IIV Bl 

B. Virtual tunneling processes (X2) 

The interface exchange term XI is a property of the separate interfaces and they contribute independently. The 
second type of exchange (X2) felt by the spins on the island is a property of the entire device. It originates from 
virtual tunneling processes, corresponding to single-electron transfer from and to the cluster. In the tunneling regime, 
this process can be treated and understood in terms of perturbation theory. In the absence of tunneling, the number 
of electrons on the island is a good quantum number. The perturbation by the contact to the electrodes allows mixing 
in of states in which the number of electrons on the island is changed by unity, at the cost of the charging energy. In 
second order perturbation theory this results in an energy gain represented by a sum over (virtually) excited states 
in which the Coulomb energy appears in the denominator and the tunneling probability in the numerator. When 
the leads are non-magnetic, these virtual processes correspond to a quantum correction to the average charge on the 
central electrode4242i This effect depends strongly on the applied gate voltage. When the unperturbed N + l (N — 1) 
particle ground state is tuned in energy just above the N particle state, the quantum correction will be large and 
positive (negative). At the degeneracy point, perturbation theory breaks down, but the ensuing divergence can be 
controlled by taking into account finite temperatures. 

When the tunneling probabilities to the ferromagnetic contacts are spin dependent, the deviations from the exact 
quantized charge on the island become spin-dependent, and therefore lead to a net excess of spins in the ground state 
that depends on the configuration of the contact magnetizations. For a symmetric spin valve it is easy to see that 
the island ground state magnetization due to these virtual processes X2 is maximal for parallel magnetizations and 
vanishes for antiparallel ones. 

The additional exchange affects non-equilibrium electron transport, in contrast to higher-order so-called co-tunneling 
processes, to the same order as the in and out-tunneling processes. For a quantum-dot island with a single quantized 
level, Konig and Martinekii showed that in the case of non-collinear magnetizations the non-equilibrium spins on an 
island injected by a finite source-drain voltage are dephased by precessing around the effective exchange field. This 
effect was also discussed for few-level quantum dots. 22 Since XI discussed in Sec. IIII Al is a material constant, the 
gate voltage dependence of X2 provides a handle for an experimental discrimination of the two effects. We derive an 
expression for the effective X2 exchange field for a classical SV-SET in Sec. IIVBI 

IV. CHARGE AND SPIN TRANSPORT 

We compute the transport characteristics of the SV-SET in lowest-order perturbation theory 44 for a diffusive or 
chaotic island in the sequential tunneling regime. The rate equations lead to a probability distribution for the excess 
number of charges n^- The excess spin accumulation s contains a large number of spins and we are interested in its 
average value in the steady state that is found from the condition (ds/dt) = 0. 

A. Charge transfer 

The operators for the excess number of electrons on the island and on the two leads are — J2ks c ks c ks 
and npa = ^2ks a aks a aks, respectively. The unpaired spin angular momentum on the cluster is written as 
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s = {H/2) 53fess' c \ s ^ss'Cks' ! where it = (er^, <7 a , cr^) is the vector of Pauli spin matrices. It is convenient to intro- 
duce a vector chemical potential Ap in the island, with size Ap — 2 |(s)| / {p^K) (see also Ref. Il3|) . where /9/v is the 

density of states at the Fermi energy. We can take into account Stoner enhancement intra-island exchange effects in 
terms of the static susceptibility \s, and we may also write Ap = 2p? B \{s)\ / {xsK)- We denote the unit vector in the 
direction of the spin accumulation by s. 

The charge current is equal to the expectation values for the rate of change of n^. In terms of the tunneling 
Hamiltonian H F = Ht\ + Hti the time-evolution is given by 



dn 



N 



dt 



- [H T ,n N ] 



— ft T kqs' aj Lks' C aqs' + h.C. 

akqs' 



(7) 



We use the interaction representation, and write the total Hamiltonian as H = H' + Ht- To second order in Ht we 
have 



dn ]v (i) 
dt 



dt' 



dnN (t) 
Jt ' 



H T (f) 



where (..) denotes an expectation value with respect to Hamiltonian H' . The electrochemical potentials of the two 
reservoirs are p C Fi = eV/2 and p C F2 = —eV/2. It is convenient to introduce grand canonical Hamiltonians including 
the chemical potentials asi£*2i 



Kn = Hn — h 1 Ap ■ s, 
K Fa = H Fa - p cFa n Fa . 

The time dependence Ck s it) = e^ KNt Ck s e'^ KNt can be formulated in terms of the projection operators 



« T (*) = H ( J + S ■ 7 ) » 



u l (s) = ^ (/ - s ■ ~a) , 



(8) 
(9) 



(10) 
(11) 



where / is the unit matrix, by making use of the equality 



eh 2 



Cps" • 



(12) 



The leads and the island are supposed to be in thermal equilibrium, so that (c\. s , Ck' S ") o — f (e/ss') Skk'ds's" , with 
Fermi-Dirac distribution / (e) = (l + e' 96 ) , where (3 is the inverse temperature. Using the expression for the matrix 
elements 



U{6 a )u s " {s)U(9j 



1 



— (1 + s s s ■ m a 

with s',s"e {t,|} = {+,-}, 



8' S' 



the rate of change of the number of electrons on the island reads 

= ^2 (G a + s"P a G a s ■ m a ) x 



/ dn N \ 
\ dt / 



n N =m as i 
F 



„Ap 



-E m -i + E m — fi cFa + s 



■ - // 

F E m — E m +i + HcFa — S —— 



where F (e) = e (l — e ^ £ ) 1 and E m = e 2 (m — CaVa/e) 2 /2C. The relation between the up and down spin conduc- 
tances (Gil anc j Gil) anc j the tunneling coefficients is G^ s = (ne 2 /h) pnPFos |Tf | , where |T S T is the value of |T fc Q ? J 2 
at the Fermi energy averaged over all the modes. pFas is the spin-dependent density of states in ferromagnet a. 



(13) 



(14) 
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In the low-bias regime considered here we can linearize Eq. (|14fl in Ap, and n C Fa- The resulting expression for 
the rate for electron tunneling through contact a, increasing the excess number of electrons from "0" to "1", is 
denoted by r^T* 1 . The analogous rate for removing one electron when um is "1" is T^ . Explicitly, we find 



G a 



Ga 



K (V, V G , An) = —IfF (E - E{) + —=-F (E - E x ) -ficF* + -^P a s ■ 7Ti a 



A/i 



Ga „, 



rr u (^, V G , A/i) = —7fF (E, - E ) - —rrF {E 1 - E Q ) 



-fJ-cFa 



2 

A/i 



P n s ■ m, 



(15) 
(16) 



Now that we have determined the tunneling rates we can write down the master equation for electron transport in 
the orthodox model. We consider a regime in which eV <C ksT -C e 2 /2C, and restrict ourselves to a gate voltage 
range for which the excess number of electrons n N alternates between "0" and "1"(0 < CqVq < e), knowing that 
the results will periodically repeat with this period. The center of the Coulomb oscillation for transitions between 
n N = "0" and "1" electrons is at C G V G = e/2. 

The steady state on the island is characterized by a constant spin accumulation (to be determined below) and the 
probabilities po and pi that there are "0" or "1" excess electrons. We have po + pi = 1. The rate equation for the 
probabilities is 



d Pn /dt = - Pn (r" +1 + r"- 1 ) + Pn+1 r n+1 - + Pn _ 1 r"- 1 ^ 1 . 

From the condition of detailed balance, pqT ^ 1 = piT 1 ^ , we find 



P0 (V,V G ,A f i) = f(Eo-Et) 



13 f (E - E,) f (E, - E Q ) 
G\ + G2 



^ ( GallcFa - PaGa—^S ■ m c 



(17) 



(18) 



The expression for the conductance of the SV-SET as a function of the spin accumulation can now be calculated and 
reads 



G(V, V G , Afi) = -epoT^ 1 + ep 1 T 1 r° 

_ G1G2 (3(Eo — Ei) 



Gi + G 2 2sinh/3(E -Ei) 



1 - 



A^i 

2eV 



s ■ (Pimi - P2W2) 



(19) 



B. Spin accumulation 



The steady-state spin accumulation is found by setting the total rate of change of s to zero. There are several 
contributions to the dynamics of the spin accumulation: 



ds\ I ds 



dt I 



7l N —0 

ds\ 



ds 



ds\ 
~dl) 



n N =l 



liiagii 



ds 
dt 



(20) 



The first two terms are due to the tunneling processes, the remaining ones to exchange, external magnetic fields and 
spin flip. We start from 



ds i r . 



(21) 



with an expectation value that to second order in Ht reads 

t 



ds(t) 
dt 



ds{t) 
dt 
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The spin current (rate of change of the spin angular momentum) due to tunneling when m excess electrons are on 
the island reads (cf. Eq. JT^l) 



ds 
dt 



= 4^2 X! ( G « s "« + P a G a m a ) x 

riN=m as" 

- F ( -E m -l + E rn - McFa + S " + F (- E m+1 + E m + HcFa ~ " 

2_\. PaG a s" (m a x s) x 

as" 

/(ei)(l-/(e 2 )) 



2 



de-i / efe 



dei / de 



e 2 - t\ + E m -i - E m + fi cFa - s"-£ 
f(e 2 )(l-f(e 1 )) 



e 2 — e i + E m — E m+1 + n c p, 



(22) 



where the prime denotes a principal value integral. Here we used the relation: 

U (9 a ) u s " (s) a l U (0 a ) f ] = i S "s + h'Jn a + -is's" (nt a x s) , 



(23) 



with s',s"e {t,|} = {+,-}. 

To first order in the small induced energy shifts the exchange Hamiltonian H cxa modifies the unpaired spins as 



ds 
dt 



[H exa ,s(t)} 



which results in a precession: 



ds 

~dt I HM 



(24) 



(25) 



where M is the number of transport channels in the normal metal and the energy shifts Ae are found in Eq. (|A2|1 . 
The conductance parameters of an F|N contact areSI 



= 6 T Y. i*™ - C (C)') - («, 6t, 4). 



(26) 



Here n and m denote the transport channels in the normal metal and r*™ and r"™ are the corresponding spin- 
dependent reflection coefficients. The contact conductances for spin-up and spin-down electrons are GjJ and G\} and 
the mixing conductance G]} governs the transverse spin currents that are absorbed and reflected by the ferromagnet 
a. The current polarized normal to the magnetization but in the plane of s and fh a is proportional to Re G]} and 
describes the spin-transfer to the magnet, thereby dissipating the spin accumulation. In the case of tunnel junctions 
ReGJ,J- — i G a /2. The out-of-the s,m a plane spin current is caused by reflection processes that make spins precess 
around m a and is proportional to ImG^. This mixing conductance has been evaluated from first principles for 
various contact materials and is small for intermetallic interfaces because positive and negative contributions in the 
space spanned by the transport channels average outM> However, there is no general reason that ImG^ should be 
smaller than G or ReG^. It is known to be quite large for the Fe|InAs interface^ and found to be very significant for 
the magnetization dynamics of MgO magnetic tunnel junctions. 2 ^ For a simple model barrier discussed in Appendix 
IbI we find the value ImG^/G = —0.26. Using the relation between the reflection phases and the energy shifts as 
derived in Eq. l|A2fl . we can rewrite the contribution given in Eq. (|25|l in terms of the imaginary part of the mixing 
conductance as (cf. Ref. |23|) 



ds\ 1-m.Gl 1 . 



dt / e xa PNC 



-m. 



(s) . (27) 



8 



The spin accumulation can also be affected by a magnetic field B, which can either be externally applied, a stray 
field from the ferromagnets, or an internal anisotropy field. The spin accumulation induced by this magnetic field can 
safely be neglected, but the induced precession of the spin accumulation is relevant, and is given by 



ds\ 
dt) 



9HB 



(28) 



magn 



Finally, spin-flip relaxation in the normal metal is taken into account by spin-accumulation decay with a spin-flip 
relaxation time r s f , 



Tat ' 



(29) 



Combining the terms in Eq. H2()(l . the spin accumulation should fulfil the stationary state condition: 

dt V ' G 7 2e 2 2sinh/?(£; -£; 1 ) 

Gi + G 2 eV {PlThl ~ P2 ™ 2) ~ (Gl + ° 2) ^( S+ ( S '^) b 

- — — B cS x (s) = 0, 

n r s f 



where 



6= m | P2G2 

Gi + G2 G\ + G; 



■m 2 . 



(30) 



(31) 



The total effective magnetic field B c g consists of the external magnetic field and contributions from the exchange 
effects XI and X2, and reads 



B cS (V G ) =B + B X1 + B X2 (V G ) , 



(32) 



Here we introduced 17 



with Bxi 

Bx2 (V G ) 



p N gfiBe- 

h 



y^Im G]}m a , 

a 

— (G 1 + G 2 )b 



i/ (E a ~ Eh.) J def (e) r, (e + E - E_ u 
- 1/ ( El - E ) j def (e) v (e + E x - E , & - 



T] (s, U) = J dm 

= - Re 



i-/M f( w ) 



m — e 



TD — £ — U 



2 2tt J V 2 2tt 



(33) 
(34) 



(35) 



where \1/ (z) is the Digamma function. In appendix [C] we discuss the derivation of the expression for B X2 i n more 
detail and comment on the differences compared to the case of a single- level quantum dot. 



V. RESULTS AND DISCUSSION 



Fig. 13 shows the magnitude of the total effective magnetic field B c g as a function of gate voltage (solid line) for a 
symmetric spin valve with parallel magnetizations (it vanishes for the antiparallel configuration) and a polarization 
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0.5 



C G V G /e 

FIG. 2: The effective magnetic field strength |-B c ff | as a function of gate voltage (solid line) for a spin valve in the parallel 

configuration. The parameters are Gi = G2, Pi = P2 = 0.7, ImG\ l = lmG\ l = Gi/4 and e 2 / (2G) = lOfcsT. The imaginary 
part of the mixing conductance gives a constant offset (dotted). The dot-dashed/dashed curves show the effective field for zero 
and one excess electron on the island. 



Pi = 0.7. The XI term is a constant that does not depend on gate voltage (dotted line). Bx2 vanishes when CgVg 
equals 0, e/2 and e. At these points, contributions from incoming and outgoing electrons cancel each other (see 
appendix[nj. The curve repeats as a function of gate voltage with period c/Cq- The spin accumulation on the island 
found from Eq. (|30[1 tends to suppress the current through the system. Spin-flip and exchange effects that dissipate or 
dephase the spin-accumulation therefore increase the conductance. As a reference we list here the conductance G (9) 
for a spin valve without interaction, with equal conductance parameters for the left and the right tunneling barrier 

Gi = G'2, P\ = P'2- 



G(0) = % 



^ PjGj (gi + 2G sf ) sin 2 9/2 ^ 



V 



Gi +2G sf + (2cosf ImG^) y 



(36) 



The final result for the symmetric spin-valve with interaction can be obtained simply from this expression by the 
substitutions: 

13 {E - E-l) 

Gl 2esinh/3(£; -P 1 ) Gl ' (37) 
n 2 cos 5 

For nonmagnetic contacts (Pi = 0) this result reduces to the known expression for normal metal single-electron 
transistors^ 

As shown in Fig. |3Ja), changing the relative strengths of XI and X2, or, since the X2 contribution is proportional to 
the polarization of the leads, Im G]} /P a G a , qualitatively modifies the current profile of the Coulomb oscillations. The 
constant offset given by Bxi skews the exchange field around CgVg = e/2, causing asymmetric conductance curves. 
When the offset starts to dominate the symmetry gets restored. The X2 contribution vanishes when the Coulomb 
blockade is lifted [CgVg = e/2), so the angular dependence of the conductance for different values of ImG^ /G\ in 
Fig. Bib) reflects only the XI effect. The curve is a simple cosine for ImG^ = 0, but is sharpened for larger ImGj^ 
because of the dephasing of the spin accumulation occurring for noncollinear angles. In Fig. [21(c) ImG^/Gi is fixed 
to 0.25 and curves are plotted for different values of the gate voltage. It can be seen that the angular dependence 
differs because the X2 depends on Vq in an asymmetric way around CqVg = e/2. 

As can be seen in Fig. 0] the shape of the Coulomb oscillation can develop minima when the polarization is high 
and the magnetizations are nearly antiparallel. At the values of gate voltage where the XI and X2 exchange effects 
cancel, the spin accumulation is not dephased and the conductance is suppressed. 

Fig. 0shows results for the conductance and spin accumulation as a function of applied magnetic field in the x (solid 
line), y (dashed) and z (dotted) directions. The spin valve is again symmetric with P\ = 0.7 and ImG}^ = G%/4. 
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FIG. 3: (a) Coulomb oscillations at fixed angle 8 — it/2 for a symmetric SV-SET with ratio lmG\ /Gi = (solid), 0.25 
(dashed) and 1 (dotted) in units of G\. The polarization P is 0.7. (b) Conductance as a function of the angle for the same 
parameters as in a), with CgVg fixed to 0.5. (c) Conductance as a function of 9 with ImG^/Gi = 0.25. Results are shown 
for CgVg equal to 0.45 (dashed), 0.5 (solid) and 0.55 (dotted). 



The angle 9 is fixed to ir/2 and CqVg — e/2. The conductance then depends only on the x-component of the spin 
accumulation (see Eq. I|19[) ). Without applied magnetic field, the spin accumulation has components in the x and 
y directions, while B c ff is in the y-direction. The results can be understood in terms of the dephasing of the spin 
accumulation by the magnetic-field induced precession that, for sufficiently large and non-collinear magnetic fields, 
quenches the spin accumulation. This "Hanle" effect is responsible for the conductance minimum at negative applied 
magnetic field in the y direction. In(5Jc) only two curves are visible because the curves for magnetic fields in the x 
and y direction overlap. 



VI. SUMMARY 



We studied the transport properties of single-electron spin valve transistors as a function of the magnetization 
configurations in the orthodox model of the Coulomb blockade. Two types of exchange effects between the spin 
accumulation on the island and the lead magnetizations play a role: a nonlocal interface exchange effect (XI) and 
exchange due to virtual tunneling processes (X2). For metallic dots these two effects are found to be of comparable 
magnitude. We predict that a line shape analysis of the Coulomb oscillation peaks should help to experimentally 
disentangle the two contributions. Additional information can be obtained by the Hanle effect. 
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FIG. 4: Conductance as a function of gate voltage for a symmetric SV-SET with 6 — 0.9-7T and ImGj /Gi = 0.15. Curves are 
shown for Pi = 0.7 (solid), Pi = 0.85 (dashed) Pi = 1 (dotted). 
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Let us consider a normal metal island in contact to a ferromagnet by a tunnel barrier (See Fig. without Coulomb 
interaction. We wish to calculate the spin-dependent shifts of the energy levels due to the presence of the F|N contact. 
In Ref. an analogous calculation was done for a ballistic one-dimensional quantum wire. Here we consider an island 
in the quasi-classical regime, i.e. the de Broglie wavelength is much smaller than the size of the island. 

The Bohr-Sommcrfcld quantization rule^i 



can be used to find the energy shifts, where p m (x) is the classical momentum for an electron in mode m, and n is 
an integer. The integral is over a whole period of the classical motion in the quasi-classical region. The total phase 
shift due to the reflections at the turning points is ifi™ + </>™, where 0™ is the spin- independent phase shift picked 
up during the reflections from the boundaries for an isolated island without contact to the ferromagnet. The small 
spin-dependent phase shift </>™ <C 1 arises from the weak coupling to the ferromagnet. The phase shifts have to be 
computed quantum mechanically via the spin-dependent reflection coefficients r™" 1 for mode m at an interface that 
is assumed to be specular (see also Appendix IB! . 

From Eq. (|A1|) . we see that increasing the quantum number n by one corresponds to introducing an extra phase 
period that increases the kinetic energy by M/ppj, where pm is the density of states of the island and M is the number 
of modes. The energy shift for an electron in mode m is therefore, to linear order in </>™, 



VII. ACKNOWLEDGEMENTS 



APPENDIX A: ENERGY SHIFTS 




(Al) 



Ac 




(A2) 



The effect of the interface on the island states can be taken into account by introducing an effective Hamiltonian as in 
Eq. J(jJ). In the case of a spin-independent tunneling barrier to a ferromagnet, the spin-splitting of the energy levels 
is small, of the same order as the transmission probability (see App. IBl . 
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FIG. 5: (a) Conductance as a function of a magnetic field applied along the the x (solid), y (dashed) or z (dotted) direction 
in units of G\. The SV-SET has symmetric junction parameters, with polarizations Pi = 0.7 and Im G[' = Gi/4. The 
magnetizations are fixed to mi/2 = (±1, 0, l)/\/2, yielding an angle 8 = ir/2 and CgVg = e/2. (b,c,d) The x, y and z 
components of the spin accumulation for the same parameters. The curves in (c) for magnetic fields in the x and y directions 
overlap. 



N F 



FIG. 6: A normal metal island with tunnel contact to a ferromagnetic lead. 



APPENDIX B: RECTANGULAR BARRIERS 

Here we evaluate the spin-mixing conductance for a model barrier, giving more details of the results of Ref. l23l 
We consider a smooth rectangular barrier between a normal metal and a Stoner- model ferromagnet. The solution of 
the Schrodinger equation for spin s in the normal metal, tp™ (x, y, z) can be used to determine the reflection coefficients 
r mm £ or gach mo de m. It reads 

V s m (x, y, z) = ^ {X lh Z) U h » x + rT m e~^) , 

V k N 

where x™ (y, z) is the transverse wave function and k™ is the longitudinal wave number for mode m in the normal 
metal. In terms of the wave numbers in the normal metal k™, barrier k 7 ^ and ferromagnet k™ s for a given energy, the 
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reflection coefficient for mode m at the barrier reads 

C m = P(kw, k%) + e"»TB, k%)p(k%, (Bl) 
where a is the barrier thickness and 

2VM2 , no x 
r(fci,fc 2 ) = (B2) 
fci + fc 2 

For a tunneling barrier, is imaginary and the spin dependent correction to the reflection coefficient is exponentially 
small in the barrier thickness. 

For a numerical estimate we use a Fermi energy in the normal metal of 2.6 eV, a barrier height of 3 eV and barrier 
thickness of a = lnm. The Fermi momenta in the ferromagnet are taken to be &ft = 1.09^4 _1 and fc^x = 0.42^4 _1 
(characteristic for Fe, see Ref. I39I) . 

For the spin-mixing conductance G^, Eq. (|26[) . we find that ImG^/G = —0.26 for this choice of parameters. The 
effective field due the interface exchange effect is therefore not negligible compared to the conductance parameters. 
More realistic electronic structure calculations^ 5 should be carried out to obtain better estimates. 

APPENDIX C: X2 EXCHANGE IN CLASSICAL DOTS 

Here we present more details concerning the derivation of Eq. (|34|l for the effective exchange field X2 in classical 
SV-SET's, that complement the derivation in Refs. ^3 f° r single-level quantum dots. Since the model is periodic 
in the gate voltage with period e/Cc, we restrict our discussion to the range < GqVq < e. From Eq. the 
contributions from virtual tunneling processes to the rate of change of s then read: 



\dt/ X2 \dt/ X2nN=0 \dt/ X2nN=1 
Using the spin currents from Eq. 1221) . we obtain, e.g.: 



ds\ h 
It 



/ = 4^ ^2 p a G a s"(m a x s) x 

/ _X2.niv— as" 

dei / de 



X2,n N =0 

/(ei)(l-/(e 2 )) 




62 — €\ + E-i — Eq + fl c Fa 

.UJ^ WHM — ^|. (C1) 

J J ( e 2 - ei + E - Ei + \i cFa - 

The first term in brackets describes virtual processes in which an electron tunnels out of the island, and the second 
term corresponds to incoming electrons. The expressions for the energy differences are given by 

£_! - E Q = (C G V G + e/2) e/C, (C2) 
Eo-E! = {CqVq - e/2) e/C. (C3) 



Because of the periodicity in the gate voltage 

/ds\ _/^£\ 



(-jt) = { -jr > with Vq > Vq — e/C G . 

\ at I x%,n N —\ \ at / X2,n N =Q 

We can now rewrite Eq. (|C1|) in terms of the function 77 (e, U), defined in Eq. I|35l) . The probabilities po and pi are 
taken from Eq. i|18|) . After linearization in V and A/i, we arrive at the expression Eq. (|34|l . 

We note the differences with the results for single-level quantum dots>i£ our expression includes an additional 
integral over the island states. For a single-level quantum dot, X2 is active only when exactly one electron resides on 
the dot, since there is no unpaired spin in an empty or doubly occupied dot. In contrast, a net spin accumulation 
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can reside on our classical dot for all numbers of electrons. The effective magnetic field is a sum weighted with the 
probabilities for "0" and "1" electrons on the dot, which leads to a partial cancellation of the contributions for different 
ri/v, as is shown in Fig. 
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